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Abstract 

The XXZ Gaudin model with generic integerable boundaries specified by generic 
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1 Introduction 



Gaudin type models £Q have played an essential role in the study of quantum systems 
appeared in many physical fields such as the BCS theory [2] of small metallic grains [HUH IE], 
nuclear physics theory Ej and QCD theory (HI E] ■ They also provide a testing ground for 
ideas such as the functional Bethe ansatz and general procedure of separation of variables 
[Tn| ITU IT2] and the construction of integral representations of the solutions to the Knizhnik- 
Zamolodchikov (KZ) equation [13 d EES EEE1 EEZ| ■ 

The original Gaudin's magnet Hamiltonians (or Gaudin operators) can be rewritten in 
terms of the quasi-classical expansion of the transfer matrices (row-to-row transfer matrices) 
of inhomogeneous spin chains with periodic boundary condition |18j . Since the elegant work 
of Sklyanin [12], the powerful Quantum Inverse Scattering Method (QISM) has been applied 
to various integrable models with non-trivial boundary conditions, which are specified by 
K- matrices satisfying the reflection equation and its dual [20J. The quasi-classical expansion 
of the corresponding transfer matrices (double-row transfer matrices) produces generalized 
Gaudin Hamiltonians with boundaries fUJ [!H] • in particular, twisted boundary conditions 
and open boundary conditions associated with diagonal K-matrices give rise to Gaudin 
magnets in non-uniform local magnetic fields ^H] and interacting electron pairs with certain 
non- uniform long-range coupling strengths [H E21 ESI EH IS] ■ 

In this paper, we study the XXZ type Gaudin magnets with most generic boundary 
conditions specified by the generic non-diagonal K-matrices given in [231 EZSj . In section 
3, we construct the generalized Gaudin operators associated with the generic boundary 
K-matrices. The commutativity of these operators follows from the standard procedure 
[T%1 ITdI |2T| specializing to the inhomogeneous spin-| XXZ open chain, thus ensuring the 
integrability of the Gaudin magnets. In section 4, we diagonalize the Gaudin operators 
simultaneously by means of the algebraic Bethe ansatz method. This constitutes the main 
new result in this paper. The diagonalization is achieved by means of the technique of 
the "vertex-face" transformation. In section 5, we conclude this paper by offering some 
discussions. 
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2 Preliminaries: the inhomogeneous spin-| XXZ open 
chain 

Throughout, V denotes a two-dimensional linear space and a , a z are the usual Pauli ma- 
trices which realize the spin-^ representation of the Lie algebra sl(2) on V. The spin-| XXZ 
chain can be constructed from the well-known six- vertex model R-matrix R(u) G End(V (E) V) 
[2*7] given by 



R(u) 



The coefficient functions read 



b(u) c(u) 
c{u) b{u) 



(2.i; 





:«) 


— -^22 


:«) 


= a(u 


R 12 




= #21 


:«) 


= &(«; 


R 21 




= #12 







1, 

smit 
sin(w + if) 
sin 77 



(2.2) 



sin(« + 77) 

Here if is the spectrum parameter and 77 is the so-called crossing parameter. The R-matrix 
satisfies the quantum Yang-Baxter equation (QYBE), 

#12(1*1 - u 2 )R 13 (u 1 - u 3 )R 23 {u 2 - u 3 ) = R 23 {u 2 - u 3 )R 13 (ui - u 3 )R 12 {ui - u 2 ), (2.3) 

and the properties, 

Unitarity : R 12 {u)R 21 {-u) = id, (2.4) 

Crossing-unitarity : (R) 2 \(-u - 2n)(R)%(u) = smusm ( u + 27 l) id (2 .5) 

sm(i* + 77) sm(u + fj) 

Quasi-classical property : R^u)^^ — id. (2.6) 

Here R 2 \{u) = Pi2#i2(w)Pi2 with P\ 2 being the usual permutation operator and U denotes 
the transposition in the z-th space. Here and below we adopt the standard notations: for 
any matrix A G End(V^) , Aj is an embedding operator in the tensor space V (g> V ® • • •, 
which acts as A on the j-th space and as identity on the other factor spaces; Rij(u) is an 
embedding operator of R-matrix in the tensor space, which acts as identity on the factor 
spaces except for the i-th and j-th. ones. 
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One introduces the "row-to-row" monodromy matrix T(u), which is an 2 x 2 matrix with 
elements being operators acting on V® N , where N = 2M (M being a positive integer), 

T (u) = Rqx(u + z x )Rq 2 {u + z 2 ) ■ ■ ■ Ron(u + z N ). (2.7) 

Here {zj\j = 1, • • • , N} are arbitrary free complex parameters which are usually called inho- 
mogeneous parameters. With the help of the QYBE ()2.3j) . one can show that T(u) satisfies 
the so-called "RLL" relation 

R 12 {u - u)r a («)T 2 (i;) = T^T^R^u - v). (2.8) 

Integrable open chain can be constructed as follows ^H]. Let us introduce a pair of 
K- matrices K~(u) and K + (u). The former satisfies the reflection equation (RE) [20] 

i?i 2 (wi - u 2 )Ki(ui)R 2 i(u 1 + u 2 )K2(u 2 ) 

= K^{u 2 )R 12 { Ul + u 2 )K{ K)7? 2 i(mi - u 2 ), (2.9) 

and the latter satisfies the dual RE 

R 12 {u 2 - u^K+^R^-ui -u 2 - 2r])K+(u 2 ) 

= K+(u 2 )R 12 (- Ul -u 2 - 2i 1 )K+{u 1 )R 2l (u 2 - Ul ). (2.10) 

For open spin-chains, instead of the standard "row-to-row" monodromy matrix T(u) ()2.7|) . 
one needs to introduce the "double-row" monodromy matrix T(u) 

T(u) = T(u)K-(u)T-\-u). (2.11) 

Using (|2.8|) and ()2.9|) . one can prove that T(u) satisfies 

R X2 {ui - u 2 )Ti(ui)R 21 (ui + u 2 )T 2 (u 2 ) = T 2 {u 2 )R 12 {u 1 + u 2 )T 1 {u 1 )R 2 x{ux - u 2 ). (2.12) 

Then the double-row transfer matrix of the inhomogeneous spin-| XXZ chain with open 
boundary is given by 

t(u) = tr(K + (u)T(u)). (2.13) 

With the help of (|2.3J) - ()2.5JI and ()2.9|) - ()2.10j) . one can prove that the transfer matrices with 
different spectral parameters commute with each other [T9*] : 

[r(u),r(«)] = 0. (2.14) 

This ensures the integrability of the inhomogeneous spin-| XXZ chain with open boundary. 
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3 XXZ Gaudin models with generic boundaries 



In this paper, we will consider a generic K-matrix K (u) which is a generic solution to the 
RE ()2.9)1 associated the six- vertex model R-matrix [2H1 EEj 

K-(u)=( %M JiS 1 = K(u). (3.1) 
The coefficient functions are 

k\{u) 

k l 2 (u) 

k\{u) 



kl{u) kl{u) 

2 cos(Ai - A 2 ) - cos(Ai + A 2 + 2£)e~ 2iu 
4 sin(Ai + f + u) sin(A 2 + f + u) 
-isin(2tt)e- i(Al+A2) e-™ 
2 sin(Ai + | + u) sin(A 2 + £ + u) ' 

2sin(2w)e i(Al+A2) e- iu 
2 sin(Ai + I + u) sin(A 2 + f + u) ' 



2 _ 2 cos(Ai - A 2 )e~ 2Mt - cos(A 1 + A 2 + 2fl 
fe2W " 4sin(A 1 + e + n)sin(A 2 + ^ + n) " 

At the same time, we introduce the corresponding dual K-matrix K + (u) which is a generic 

solution to the dual reflection equation (j2.10j) with a particular choice of the free boundary 

parameters with respect to K~(u): 

k + (u) = ( ?M ) • (3.3) 

The matrix elements are 
k+\(u) 

k + l(u) 



u 



k+t(u) k+l(u) 

2 cos(Ai - \ 2 )e- lT > - cos(Ai + A 2 + 2()e 2lu+lT i 
4 sin(Ai + ^ — u — rj) sin(A 2 + £ — u — rj) 
i sin(2u + 2rj)e-^ Xl+X2 h iu - in 
2 sin(Ai + £ — u — rj) sin(A 2 + £ — u — rj) ' 

-i sin(2w + 2r/)e i(Al+A2) e iM+ir ? 
2 sin(Ai + £ — u — rj) sin(A 2 + £ — u — rj) ' 



a , +2(m) = 2 cos(A 1 - A 2 )e 2Mt+ "? - cos(A 1 + A 2 + 2j)e~^ 
4 sin(Ai + £, — u — rj) sin(A 2 + ^ — u — rj) 

The K- matrices depend on three free boundary parameters {Ai, A 2 , £} which specify inte- 

grable boundary conditions j2S|. It is very convenient to introduce a vector A = Ylt=i ^k^k 

associated with the boundary parameters {A.;}, where {q, i = 1,2} form the orthonormal 

basis of V such that (q, €j) = 5ij. We remark that K~{u) does not depend on the crossing 

parameter n but K + (u) does. They satisfy the following relation: 

\im{K + (u)K-(u)} = \im{K + (u)}K(u) = id. (3.5) 
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Let us introduce the generalized XXZ Gaudin operators pQ {Hj\j = 1,2, ■ ■ ■ , N} as- 
sociated with the spin-| XXZ model with generic boundaries specified by the boundary 
K- matrices in ()3.1|) and (|3.3|) : 



2M r . , 



#i = r i(^) + sin ( z .- Zk ) \ a K a i + ^ a t + cos ^' " 

2M r^— 1/ \ r z z 1 



- 1 



where T^u) = ^{^(«)}|^=o^(«), J = l,---,iV, with = tr {^ + (n)i? 0j (2n)F i}, 

and {zj} correspond to the inhomogeneous parameters of the inhomogeneous spin-| XXZ 
chain with generic open boundary. For a generic choice of the boundary parameters {Ai, A2, £}, 
Fj{u) is an non-diagonal matrix, in contrast to that of 21J. 

The XXZ Gaudin operators (J3.6)) are obtained by expanding the double-row transfer 
matrix t(u) ()2.13|) at the point u = Zj around 77 = 0: 

T(z 3 )=T(z 3 )\ v=0 + r ] H J + O(r ] 2 ), j = l,---,N, (3.7) 
H i = -^T(zi)\n=o- (3-8) 

The relations (|2.fij) and ()3.5|) imply that the first term r(zj)\ v= o in the expansion ()3.7j) is 
equal to an identity, namely, 

r(z i )U=o = id. (3.9) 

Then the commutativity of the transfer matrices {t(zj)} (j2.14|) for a generic r] implies 

[H j} H k ] = 0, z,j = l,---,N. (3.10) 

Thus the Gaudin system defined by ()3.6j) is integrable. Moreover, the fact that the Gaudin 
operators {Hj} ()3.6|) can be expressed in terms of the transfer matrix of the inhomogeneous 
spin- 1 XXZ open chain enables us to exactly diagonalize the operators by the algebraic 
Bethe ansatz method with the help of the "vertex-face" correspondence technique, as can 
be seen in the next section. The aim of this paper is to find the common eigenvectors and 
eigenvalues of the operators ()3.6|) . 
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4 Eigenvalues and Bethe ansatz equations 

4.1 Six- vertex SOS R- matrix and face- vertex correspondence 



The simple root a and fundamental weight Ax of si (2) are given in terms of the orthonormal 
basis {ej} as: a = e\ — e 2 , Ai = ~. Set 

2 2 

2 

fc=l i=l 

For each dominant weight A = aAi , a G Z + (the set of non-negative integer), there exists an 
irreducible highest weight finite-dimensional representation V\ of Ai with the highest vector 
| A). For example the fundamental vector representation is Va 1 . 

Let f) be the Cartan subalgebra of Ai and h* be its dual. A finite dimensional diago- 
nalisable f)-module is a complex finite dimensional vector space W with a weight decom- 
position W = ©^6fj*W'[/x], so that I) acts on W[fj] by xv = fJ,(x)v, {x G f), u G 
For example, the fundamental vector representation V\ 1 = V, the non-zero weight spaces 
W\i\ = Cu, i = 1,2. 

For a generic to G V, define 

m, = (m, ej), TOjj = TOj TOj = (to, ej - e^), £, j = 1, 2. (4.2) 

Let -R(m, to) G End(V^ ® V) be the R-matrix of the six-vertex SOS model, which is trigono- 
metric limit of the eight-vertex SOS model given by 

2 2 

R(u, m) = J2 Ki% m)E tt ®E ii +J2 m)Eu Ejj ■ flg(u, m)E Jl ®E ij } , (4.3) 

i=l 

where is the matrix with elements {Eij) l k = SjkSu. The coefficient functions are 

*S(«,A)=1, flg( M ,A)= l+h ( 4 . 4) 

J sin(u + r]) sm(TOy) 

civ \ sin7ysin(M + TO i3 -) 

Rij{u,m) = — — r, z ^ J, 4.5 

J sin(u + 77) sm(TOjj) 

and TOjj is defined in (14. 2 j) . The R-matrix satisfies the dynamical (modified) quantum Yang- 
Baxter equation 

Hi 2 («i - «2, to — r]h^)R 13 (ui - w 3 , m)R 2 3{u 2 — u 3 , m — rjh^ 1 ') 

= R2z{u 2 - u 3 ,m)R 13 (u 1 -u 3 ,m- T}h ( ^)Ri 2 {ui - u 2 ,m). (4.6) 



We adopt the notation: Rn(u, m — rjh^ 3 ') acts on a tensor v\ ® i>2 <8> V3 as R(u, m — r]jj) ® id 
if v 3 e W\p\. 

Define the following functions 

eV*(u) = e- iu , 6 {2 \u) = l. (4.7) 

Let us introduce an intertwiner, i.e. an 2-component column vector 4 > mm-m( u ) wnose k-th 
element is 

0« Ju) = eto(u + 2m i ). (4.8) 

m,m—r]j v ' 

Using the intertwiner, one can derive the following face-vertex correspondence relation [2H] 

m,m~rfi (^1.) ® 0m— r]i,m— (^2) 
= E - «2, m)S^-^,A- n ((+fc)( tt l) ® ^ m ,m-4^)- ( 4 - 9 ) 

Then the QYBE f!2.3|) of for the vertex-type R-matrix R(u) is equivalent to the dynamical 
Yang-Baxter equation ()4.6|) of the SOS R-matrix R(u, m). For a generic m, we can introduce 
other types of intertwiners 0, satisfying the conditions, 

2 

XI C.m-wM ^m,m-r,i>{ U ) = <W> ( 4 - 10 ) 
fe=l 

EfiUwtfUw-^. ( 4 - n ) 
fe=i 

from which one can derive the following relations: 



E ^£ m ->» = <%> ( 4 - 12 ) 

2 

E Sn,Am(«) *2U,m(«) = ( 4 - 13 ) 

Through straightforward calculations, we find the K-matrices K ± (u) given by (|3.1|) and 
()3.3|) can be expressed in terms of the intertwiners and diagonal matrices JC(X\u) and JC(\\u) 
as follows 

K-{u)t = E^Vi), wWI^a-^N. ( 4 - 14 ) 
K» f s = E ^A^-^I^J^tM,, (4-15) 



Here the two diagonal matrices JC(X\u) and K.(X\u) are given by 

rv\i ^-tv fvf\\ ^ i-ai n tv / Sin(Ai+j-n) sin(A 2 + £-^ 

/C A u) = Drag A w 1, fc A « a = Diag — — — — — — -, . , 4.16 

sm(Ai + 4 + u) sm(A 2 + £ + u) 



IC(\\u) = Diag(fc(A|u)i, k(X\u] 



2) 



_ Diag , sin(Ai 2 -r7) sin(Ai+£ +u+rj) sin(Ai 2 +r/) sin(A 2 +^+n+?7) ^ 
sin A12 sin(Ai + ^ — u — rf) ' sin A12 sin(A2 + £ — w — ??) 

Moreover, one can check that the matrices K{X\u) and JC(\\u) satisfy the SOS type reflection 
equation and its dual, respectively [HO]. Although the K-matrices K ± {u) given by ()3.1|) and 
()3.3|) are generally non-diagonal (in the vertex form), after the face- vertex transformations 
(|4.14|) and (|4.15jl . the face type counterparts JC(X\u) and K,(X\u) simultaneously become 
diagonal. This fact enables us to apply the generalized algebraic Bethe ansatz method 
developed in [3T] for SOS type integrable models to diagonalize the transfer matrices t(u) 
(EH- 

4.2 Algebraic Bethe ansatz 

Using the relations (|4.12|) and f)4.13|) . the decomposition of K + {u) ()4.15|) and the diagonal 
property (|4.17|) . the transfer matrix t(u) ()2.13|) can be recasted into the following face type 
form 

t{u) = tr(K + (u)T(u)) 

fJt,U 

= y^0A,A-^(-^)^ + (^)0A-»7(A-t>),A^A( M ) <fi\-v(fi-i>),X-Vv( U ) T (m)0A,A-w( _M ) 

= £ )C(\\ur u T(\\u); i = J2 HM u )i*'T{\\u)£. (4.18) 
Here we have introduced the face-type double-row monodromy matrix T(X\u) 

T (M U )^ = 0A-r,(A-£>),A-^( U ) T(M)0 AiA _ r//t (-M) 

= E^a-^a-.aW T(«)^_ A (-«). (4-19) 

This face-type double-row monodromy matrix can be expressed in terms of the face type 
R-matrix R(X\u) (|4.3|) and K-matrix K.(X\u) (|4.16|) (for the details, see equation (4.19) of 
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PT]). Moreover from ([2.12)1 . f)4.9|) and (j4.13|) one can derive the following exchange relations 
among T(\\u)^: 



E E R ^ - + v(h + i2)|«i)i 



U,«2 Jl J2 

x jR(wi + u 2 , \)f 2 f 3 T(X + r)(j 3 + h)\u 2 )f 3 



= E E T ( A + + 5o)l«a)£fl(«i + « 2 , A)£* 

H ,«2 jl J2 

x T(A + r/(j 2 + ? 2 )| M i)^( Wl - u 2 , X)f%. (4.20) 
As in let us introduce a set of standard notions for convenience: 1 

A(\\u) = T(\\u){, fi(A|«) = ?^2, (4.21) 

sm^A 12 j 

D(A| U ) = Sm ( A 12 + ^) {r(A|t|) 2 _ R{2u ^ x + ^1)21^1^} (4 22) 

sm^Ai2; 

Hereafter, we adopt the convention: Ay = Aj — Aj, introduced in (|4.2)l . After tedious 
calculations, we find the commutation relations among ^4(A|m), £>(A|-u) and £>(A|-u). Here we 
give the relevant ones for our purpose, 

A(X\u)B(X - 2 V 1\V) = Mu + v)Mu-v- V ) _ j _ £ 

sin(« + f + 77) sin(u — f ) 

sin 77 sin 2u sin(w — t> — Ai 2 + rj) „. , „ , 

/; --B(A-2r/l|w)^l(A-2?7l|w) 



sin(w — i>) sin(2i> + 77) sin(Ai 2 — rf) 
sin 77 sin(tt + v + A 2 i + 2t/) 



sin(w + t> + rf) sin(A 2 i + 77) 



B{\-2r]l\u)V{\-2r]l\v), (4.23) 



V{X\u)B{X - 2 v i\v) = Mu-v + v)sHu + v + 2 v ) B{x _ j _ 2r]llu) 

sm(u — v) sm(u + v + rj) 
sinr7sin(2n + 2,) sin(u - v + A 12 - V ) _ _ ^ 



+ 



sin (it — v) sin(2w + rj) sin(Ai 2 — 17) 

sin rj sin 2v sin(2w + 2rj) sin(it — v + X12) 



sin(it + v + rj) sin(2t> + rf) sin(2n + rj) sin(Ai 2 — rf) 

xB(X-2r]i\u)A(X-2T)i\v), (4.24) 

B(X - 2r)l\u)B{X - irji\v) = B(X - 2r]l\v)B{X - ir]l\u). (4.25) 

Here we have used the identity 2 = — 1 which can be derived from (|4.1jl . 

1 The scalar factors in the definitions of the operators B(\\u) and T>{X\u) are to make the relevant com- 
mutation relations as concise as l|4.23[) - (|4.25[) . 
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In order to apply the algebraic Bethe ansatz method, in addition to the relevant com- 
mutation relations (j4.23|) - (j4.25|) . one needs to construct a pseudo- vacuum state (also called 
reference state) which is the common eigenstate of the operators A, T> and is annihilated by 
the operator C. In contrast to the case of the spin-1 XXZ open chain with diagonal K^iu) 
|19j . for the open chain with generic non-diagonal K- matrices (13.1)1 and ()3.3|) . the usually 
highest-weight state 

is no longer the pseudo-vacuum state. However, after the face- vertex transformations (|4.14jl 
and ()4.15|) . the face type counterparts K-matrices JC(X\u) and fC(X\u) simultaneously become 
diagonal. This suggests that one can transfer the spin-| XXZ open chain with generic non- 
diagonal K-matrices into the corresponding SOS model with diagonal K-matrices )C(X\u) and 
JC(X\u) given by ()4.16|) - ()4.17|) (possibly after some local gauge transformations [2H])- Then 
it is easy to construct the pseudo-vacuum in the "face language" and use the generalized 
algebraic Bethe ansatz method [31J to diagonalize the transfer matrix. We shall develop this 
game in the following. 

Let us introduce the corresponding pseudo- vacuum state \Q) 

= 0A-(JV-l)r?i,A-iVr 7 i(- 2 l) ® <P\-(N-2) v i,\~(N-l) v i(- z 2) ' " " ® <t> \,\- v l(- *n) ■ (4.26) 

The state does only depend on the boundary parameters {Ai,A2} and the inhomogeneous 
parameters {Zj}. Using the technique developed in [3T], after tedious calculations, we find 
that the pseudo-vacuum state given by (|4.26jl satisfies the following equations, as required, 

AX-N^ ^ttt 1 ^ (4 ' 27) 

V(X - N^um = ^^i + Z + u + V )sin(X 2 + Z-u-ri) 

sin(2« + 77) sin(Ai + £ + u) sin(A 2 + £ + u) 

sin(-u + z k + rf) sm(u — z k + rj) J 
C(X - Nrji\u)\n) = 0, (4.29) 
B(X - Nrji\u)\Q) ^ 0. (4.30) 

Then the so-called Bethe states can be constructed by applying the creation operator B on 
the pseudo-vacuum state 

\v u ■••,%) = B(X-2rji\v 1 )B(X-Ar]i\v 2 )---B(X-2Mr]i\vM)\n). (4.31) 
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x 



From (|4.18p . (|4.21|) and (|4.22|) we can rewrite the transfer matrices r(u) (|2.13|) in terms 
of the operators A and V 

, . sin(A 2 + i - u) sin(Ai + £ + u) sin(2w + 2ij) .... . 

r W = ■ r\ i c \ • / \ -7 r— 77^ — ; — r-4( A W 

sm(A 2 + 15 — u — Tj) sm(Ai + c, — u — r]) sm(2u + 77) 

sin(A 2 + £ + + 77) 
sm(A 2 + 4-77-77) 

Acting the above expression of the transfer matrices on the Bethe states \v\, ■ ■ ■ , %) (j4.31j) 
and repeatedly using the relevant commutation relations (|4.23|) - (j4.25|) . we obtain 

sin(2« + 2rj) sin 77 



t{u)\vi, ■ ■ ■ ,v M )=t(u)\vi, ■ ■ - jVmH 



sin(A2 + £ — u — rf) sin(Ai + £ — u — 77) 



M 



x 



^ F a \vi, Vg-!, u, Vg+x, v M ) \ ■ (4.33) 



a=l 



Here, the term associated with function t(u) is the so-called wanted term which gives rise to 
the eigenvalues and the term associated with {-F a |a = 1, • • • , M} is the so-called unwanted 
term. They are given, respectively, by 

sin(A 2 + £ — u) sin(Ai + £ — u) sin(2w + 2rj) 



t(u) 



+ 



sin(A 2 + £ — u — rf) sin(Ai + £ — u — rf) sin(2w + rf) 

-i-r sin(?j + Vk) sin(w — V). — rf) 
sin(w + v k + rf) sin(w — Vk) 

sin(A 2 + £ + u + rf) sin(Ai + £ + u + rf) sin 2u 
sin(A 2 + £ + u) sin(Ai + £ + u) sin(2w + rf) 

t \ sm(u + Vk + 2rf) sin(w — Vk + rf) 

X 



n 



fc=i 

2M 

X 

h 



sin(u + v k + 77) sin (it - v k ) 
sin(u + sin(it — z k ) 



sin(it + z fc + 77) sin(it — z k + 77) ' 



and 



_ sin 2u Q sin(A 2 + £ - i> Q ) sin(Ai + £ - i> a ) 
sin(2i; a + 77) sin(it + i>q, + 77) sin(it — v a ) 

tt sin(?j a + sin(u a - i> fc - 77) 
sin(i; a + v k + rf) sm(v a - v k ) 

sin(A 2 +£+'u Q ,+77) sm(\ 2 +£—Vg—r]) sin(Xi+^+v a +rf) sm(\\+£—v a —rf) 
sin(A 2 + i + v a ) sin(A 2 + £ - v a ) sin(Ai + £ + u a ) sin(Ai + £ — v a ) 



sin(i; Q + v k + 2rf) sin(t> Q - v k + 77) 



M 

J! sin(i> a + v k + rf sin(<u Q - v k ) 
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nsm(v a + z k ) sm(v a - z k ) \ 
— — ^ , , a = l,---,M. 4.35 

k=i sm(v a + z k + 77) sm(v Q - z k + r]) j 

The relation ()4.33j) tells us that the Bethe states \v%, ■ ■ ■ , %) are eigenstates of the transfer 
matrices t(u) if the unwanted terms vanish: F a = 0, a = 1, • • • , M. This leads to the Bethe 
ansatz equations, 

sin(A 2 + £ + u«) sin(A 2 + £ - f a ) sin(Ai + £ + u Q ) sin(Ai + £ — f a ) 
sin(A 2 +£+v a +r/) sin(A 2 +£-v Q -?7) sin(Ai+£+v Q +r/) sin(Ai+£— v a — 77) 

_ -TT sinQ a + v k + 2rj) sinQ a - v k + 77) 
1£ sin(v a + sm(v a -v k -rj) 

nsm(f Q + z fc ) sm(u Q - z fc ) 
— — — - T , a = l,---,M. 4.36 
fe=i siii(d q + z fc + 77) sm(v a - z k + 77) 

We have checked that the Bethe ansatz equations indeed ensure that the eigenvalues (|4.34j) 
of transfer matrices r(u) are entire functions. Our result recovers that of for the very 
special case z k — 0, k — 1, • • ■ , 2M. 

4.3 Eigenstates and the corresponding eigenvalues 

The relation (|3.8|) between {Hj} and {t(zj)} and the fact that the first term of (|3.7|) is a 
c-number enable us to extract the eigenstates of the Gaudin operators and the corresponding 
eigenvalues from the results obtained in last subsection. 
Let us introduce the state 

l^ (0) ) = ( ! J®-"®^ 6 1 j (4.37) 

This state can be obtained from the pseudo- vacuum state \Q) (|4.26|) by taking the limit: 
|0/ )) = lim^o 1^)- Let us introduce the matrix C(u) E End(y) associated the intertwiner 
vector (f) 

/ -i(u+2X 1 ) -i(u+2A 2 ) \ 

C(u)=[ e 1 x J, (4.38) 

and the associated gauged Pauli operator o~~(u) E End(V), 

o-{u) = C{u)(J-C{u)- 1 . (4.39) 
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Then we define the states ^{v\, • • • , %) and $ a : 



tt ( f sin(Ai + £ - u Q ) sin(w Q - ^ + A i2 ) 

»(«.,- ■ ■ , «m) = n (£ ( sin(Ai+?+ .„ t>)sinK _, i) 



sin(A 2 + £ - w Q ) sin(v a + z k - A 
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sin(A 2 + £ + v a ) sin(v a + z k ) 



) l^ (0) ), (4-40) 



' sin(Ai + i - vp) sm{vp - z fc + Ai 2 ) 



*.=n e 

/3^ a \fc=l 



sin(Ai + i + vp) sin(t>g - z fe ) 
5ill(A + f „ ) y + z A 12 ) | A m 

sm( A 2 + 4+^) sm^ + z k ) ) J 



Noting the relations ([3.8)1 . ()3.7|) and (|4.33jl - (|4.35jl . and using the same method as in [To] , 
we derive the following relations 

M 



HjVivi, ••■,%) - • • • , v M ) + ^ 1 ~ 



sin 2zj sin 2i> 



sm(Zj + v a ) sm{Zj — v a ) 

sin(A 2 + £ - t> a ) sin(Ai + £ - u a ) _, 1 
sm(A 2 + £ - Zj) sm(Ai + £ + %) J J 

The functions Ej and / Q are 

2 M . 2 

B, = cot 2Z, + £ cot(A J + { - Zs ) + £ - I — ^ , (4.43) 

2 , 2M 



f =v - + V- 

^-^ sWX • -I- P — D . ^ eirifX , -I- P -I- it A ^ «i 



^ sin(Aj + £ - v a ) sin(Aj + £ + u a ) ^ sin(w a + sin(w Q - z k ) 

M 

~ 2 Y—, r-^ V (4-44) 

^ sm(v a + v k ) sm(v a - v k ) 

The equation (|4.42|l suggests that the state *&(vx,- ■■,%) is an eigenstate of the Gaudin 
operators {Hj} if {vj\j = 1, • • • , M} is set to satisfy f a — for a — 1, • • • , M. This leads to 
the corresponding Bethe ansatz equations 

2 , 2Af 



(\ . 4- £ — i) \ a\-n( \ . 4- t 4- 11 < * <si 



^ sin(A i + i - v a ) sin(A i + £ + v a ) ^ sin(v a + z k ) sin(v a - z k ) 

M 

2 E w I wl V' a = 1 , • • ■ , M. (4.45) 

^ sm(u Q + v k ) sm(v a - v k ) 
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5 Conclusions 



We have studied the XXZ Gaudin model with generic boundaries specified by the non- 
diagonal K-matrices K (u), ()3.1|) and ()3.3j) . In addition to the inhomogeneous parameters 
{zj}, the associated Gaudin operators {Hj}, ()3.6|) . have three free parameters {Ai,A2,£}, 
which give rise to three-parameter (Ai, A2, £) generalizations of those in 0] and two-parameter 
(Ai,A2) generalizations of those in ^HJ EI]- As seen from section 4, although the "vertex 
type" K-matrices K ± {u) (|3.1j) and (|3.3j) are non- diagonal, the compositions, (|4.14j) and 
(|4.15jl . lead to the diagonal "face- type" K-matrices, (|4.1fij) and (|4.17jl . after the face- vertex 
transformation. This enables us to successfuly construct the corresponding pseudo-vacuum 
state \ fl) (|4.2fi|) and apply the algebraic Bethe ansatz method developed in 31 for the SOS 
type models to diagonalize the transfer matrix r(u) of the inhomogeneous spin-| XXZ open 
chain with the generic non-diagonal K-matrices. Furthermore, we have exactly diagonalized 
the Gaudin operators {Hj}, and derived their common eigenstates ()4.4()|) and eigenvalues 
(j4.43j) as well as the associated Bethe ansatz equations (|4.45jl . 

One may choose the K- matrix K + {u) whose matrix elements are given by 

2 cos(Ai - A 2 )e- ir ' - cos(Ai + A 2 + 2()e 2iu+ir > 



k + \{u) 
k + \{u) 
k +2 Au) 



4 sin(Ai + £ — u — n) sin(A2 + £ — u — rj) 
i sin(2u + 2rj)e- i{Xl+X2) e iu ~ iv 
2 sin(Ai + £ — u — rj) sin(A2 + £ — u — rj) ' 

-isin(2w + 2rj)e i{Xl+X2) e iu+iri 
2 sin(Ai + £ — u — rj) sin(A 2 + £ — w — 77) ' 



, +2 . . 2 cos(A 1 - A 2 )e 2iU +"? - cos(X 1 + A 2 + 2£)e-*> ( . . , 

k 2\ u ) = 7 f \ 7 f \ — • (5.1) 

4 sin(Ai + £ — it — rj) sin(A2 + £ — u — rj) 

Taking £ = £, the above K- matrix recovers that given by (|3.4j) . We have used the same 
method in subsection 4.2 to diagonalize the transfer matrix t(u) (|2.13|) with the non-diagonal 
K-matrices K~~(u) (J3.2|) and K + (u) (|5.1|) . The corresponding eigenvalue t(u) is given by 

sin(A2 + £ — u) sin(Ai + £ + u) sin(Ai + £ — u) sin(2« + 2rj) 



t(u) 



+ 



sin(A2 + £ — u — 77) sin(Ai + £ — u — rj) sin(Ai + £ + u) sin(2w + rj) 

n' 1 sin(w + v k ) sin(-u — v k — rj) 
-r-) : — -, r 
f sm[u + Vk + rj) sm[u — v k ) 

sin(A2 + £ + u + rj) sin(Ai + £ + u + rj) sin(A2 + £ — u — rj) sin 2u 
sin(A 2 + £ — u — rj) sin(Ai + £ + u) sin(A 2 + £ + u) sin(2u + rj) 
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t-t sin(w + Vk + 2rj) sin(w - v k + rj) 
J- \ sin(w + Vk + rj) sm(u - v k ) 

k=l 

x TT S m (^ + gfc) S m (^ - Z k ) , h2 , 

J-jL sin(w + Zk + rj) sin(u — Zk + 77) ' 
and the associated Bethe ansatz equations are given by 

sin(A 2 + £ + f a ) sin(A 2 + £ - v Q ) sin(Ai + £ + u a ) sin(Ai + £ — v a ) 
sm(\ 2 +£+v a +rj) sm(X 2 +^-v a -rj) sin(\ 1 + £+v a +rj) sin(Ai+£— v a — rj) 

_ tt sm(v a + v k + 2rj) sin(v a - v k + rj) 
J- J- sin(?j Q + v k ) sin(v a - v k - rj) 

k^a 

nsm(f a + z k ) sm(v a - z k ) 
— — f , — -, a = l,---,M. 5.3 

k=i sm(v a + Zk + rj) sin(v Q - z k + rj) 

For the £ = £ case, ()5.2|) and ()5.3|) result in ()4.34|) and ()4.36|) respectively. 
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